Energy dependent wavelength of the ion induced nanoscale ripple 
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Wavelength variation of ion beam induced nanoscale ripple structure has received much attention 
recently due to its possible application in nanotechnology. We present here results of Ar + bombarded 
Si in the energy range 50 to 140 keV to demonstrate that with beam scanning the ripple wavelength 
increases with ion energy and decreases with energy for irradiation without ion beam scanning. 
An expression for the energy dependence of ripple wavelength is proposed taking into simultaneous 
effect of thermally activated surface diffusion and ion induced effective surface diffusion. 

PACS numbers: 68.35.Bs, 79.20.Rf, 64.60.Ht 
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Formation of periodic undulations or ripple like fea- 
tures on various materials with typical wavelength rang- 
ing from about 10 nm to 1/im, obtained by obliquely 
incident ion bombardment, has become an active re- 
search subject due to its possible technological applica- 
tions, as varied as optical devices, templates for liquid 
crystal orientation and strain-free patterned substrates 
for heteroepitaxial growth of quantum wires. It is also 
expected that systematic study of ion beam induced nano 
ripple formation will help us to understand the basic pro- 
cesses prevalent in formation of sand dune like structures 
in deserts. Although this ion induced phenomenon was 
reported first in 1960s and then in 1970s @, §], the 
improvement in experimental conditions such as, better 
vacuum and ion beam parameters and improved surface 
characterizing tools, has enabled us to control the growth 
of these ripple like features §, |, | @, |, §, [u], The 
first widely accepted theoretical approach describing the 
process of ripple formation due to ion bombardment was 
developed by Bradley and Harper (BH) jl^]. This linear 
theory fl^| predicts the ripple wavelength and orienta- 
tion in agreement with numerous experimental studies. 
However, this theory cannot explain a number of experi- 
mental observations, such as the saturation of the ripple 
amplitude |)| , the appearance of rotated ripples |llj and 
kinetic roughening E^). Moreover, according to the BH 
theory ripple wavelength should decrease with ion energy 
but this prediction has not been confirmed experimen- 
tally so far 

Recently a formalism |l4|, |l5[ [T(| based on nonlin- 
ear continuum theory has been developed to understand 
these experimental observations not predicted by linear 
theory. In this new formalism, not only nonlinear and 
noise terms were included in the equation of height pro- 
file for eroded surface but also existence of two differ- 
ent surface diffusion processes were recognised. Based 
on Sigmund's theory of sputtering [0, the height evolu- 
tion h(x,y,t) of an ion eroded surface according to this 
nonlinear theory jl4|, |l5|, [l6| can be described by 

d t h = vV 2 h - (D T + D 7 )V 4 /i + ^(V/i) 2 + y, t). (1) 
Here v is the roughening prefactor (also called effec- 



tive surface tension) , generated by the ion bombardment 
sputter-erosion process. D T & D 1 are designated as ther- 
mally activated surface diffusion constant and ion in- 
duced effective surface diffusion (ESD) constant. The 
nonlinear term A represent the slope dependent erosion 
rate, where (V/i) define the local slopes and rj(x, y, t) 
is an uncorrelated white noise with zero average, mim- 
icking the local random fluctuation of the incident ion 
flux. All the coefficients of Eq. (1) can be made direc- 
tional in (x, y) space and their detailed expressions are 
known |l4|, [ll|. It is interesting to note that beyond a 
cross-over time, r ~ (D/v 2 ) ln(i//A), the nonlinear terms 
become dominating and this effect is expected to gener- 
ate many novel phenomena like ripple rotation in oblique 
ion incidence fllf| and formation of dots and holes for 
normally incident ion beams [fbsf . On the otherhand for 
time, f<r, nonlinear and noise terms can be neglected 
(A = 0, ?7 = 0) and Eq. 1 reduces to the BH theory 
|l2| except for the fact that surface diffusion constant K 
has two components K = D 1 + D T . In the linear regime 
the balance of the unstable negative sputter-erosion term 
(— \v\d 2 h), trying to roughen the surface and the positive 
surface diffusion term (KV h) acting to smooth the sur- 
face, gives rise to observable ripple with wavelength 



U = 2ny/2K/\iy l \ 



(2) 



where, i refers to the direction x or y along which the 
associated \vi\ (\v x \ or \u y \) is largest. 

The ion induced ESD process predicts increase of rip- 
ple wavelength (I) with ion energy (E) but the thermal 
diffusion process, assumed to be dominating at high tem- 
perature, predicts reduction oil with increasing E. While 
the increase of / with E has been observed experimentally 
H 0, [l8] , direct observation of decreasing I with increas- 
ing E has not been reported so far to the best of our 
knowledge, even when the sample temperature is varied 
|| . We have performed Ar + bombardment on Si(001) 
in 50-140 keV energy region and kept the bombardment 
process well within the linear region to verify experimen- 
tally the existence of these two diffusion processes and to 
understand the interplay between them in determining 
the wavelength of nanoscale ripple structure. 
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The ion bombardment of Si targets (p type Si(OOl) sin- 
gle crystal wafer) was undertaken with a high current ion 
accelerator (Danfysik) which is described in detail else- 
where These Si samples were irradiated in the im- 
plantation chamber (having a base vacuum in the range 
of 10~ 7 mbar) by a focused (with typical beam spot size 
of 1.5 mm and 2.5 mm) 40 Ar + beam at 60° angle of ion 
incidence with respect to surface normal of the samples. 
With a magnetic x — y sweeping system the focused beam 
was scanned over the sample to maintain homogeneous 
irradiation. The scanned area could be varied from 5 
mm x 5 mm to 15 mm x 15 mm. To check the influence 
of the beam scanning, several samples were irradiated 
without beam scanning. The beam flux in the present 
experiment was around 140 /xA/cm 2 for irradiation with- 
out beam scanning. As the area of the bombarded spot 
on the sample for unscanned focused beam was smaller 
than that for the case of scanned beam, the period of 
bombardment was adjusted to give same dose of 10 18 
ions/cm 2 for all the irradiated samples. The maximum 
period of bombardment in the present case was always 
less than 1 hr, which is far less than the estimated cross- 
over time r (about 3 hrs for 50 keV and 4 hrs for 100 
keV) |t5| ensuring the present experimental condition to 
be well within the linear regime. Thus, a higher dose of 
(3-4) xlO 18 ions/cm 2 is required to expect any nonlinear 
effect whose signature has indeed been observed in our 
preliminary results reported elsewhere ]2C| for the case 
of morphology developed on GaAs bombarded by 60 keV 
Ar beam. 

After irradiation, the samples were investigated by 
atomic force microscopy (AFM) in contact mode un- 
der ambient condition as described earlier |l9|, [2(], pl| . 
Fig. 1 shows representative AFM images of Si(001) sam- 
ples bombarded at 50 and 100 keV with and without 
beam scanning. Ripple topographies are clearly visible 
in all the AFM micrographs and they are oriented per- 
pendicular to the ion-beam projection (indicated by ar- 
rowmarks) onto the surface as predicted by linear theory. 
The ripples grew uniformly on the whole bombarded spot 
for irradiation with beam scanning whereas ripples were 
developed at the central region of the bombarded spot 
for irradiation without scanning. The mean amplitude 
of the ripples with and without beam scanning, obtained 
from AFM measurements, are also shown in fig. 1. The 
most probable observed wavelength (I) of ripples, ob- 
tained from AFM images and Fourier analysis is given in 
fig. 2. With beam scanning I grows from about 700 nm to 
about 1000 nm as the incident ion energy (E) increases 
from 50 keV to 140 keV. However, the / decays from 600 
nm to about 400 nm with the increasing E when the rip- 
ple is formed without beam scanning. If thermally acti- 
vated surface diffusion is the dominant process for surface 
smoothing, the wavelength dependence on energy (E) of 
bombardment is given by@ (~ E- 1 / 2 . On the oth- 
erhand in the absence of thermal diffusion, ion-induced 



ESD becomes dominating contributor in K and we ex- 
pect [[l4| the energy dependence as I ~ E. The variation 
of I with E of the present experiment has exhibited both 
the trends : one increasing and other decreasing with the 
energy giving rise to exponents of 0.45 and -0.56 respec- 
tively (refer Fig. 2). 

When the sample is bombarded with a focused beam 
without scanning, the beam intensity at the center of 
the bombarded spot becomes more than its peripheral 
region causing a local temperature rise. As the semicon- 
ductor sample in the present case is clamped on a copper 
block, good thermal contact is not ensured favouring the 
local temperature rise due to ion beam heating. Under 
this situation, the rise of temperature, T, from the initial 
ambient temperature , T Q = 25°C, at time t(sec) of the 
commencement of the bombardment can be calculated 
from the formula 0] 

t = f T SK f dT, (3) 

J To P - (T 4 — T 4 ) V ' 

where S,k and p are specific heat (0.836 J/g/°C), thick- 
ness (0.47 mm) and density (2.32 g/cm 3 ) respectively 
of the Si wafer used in the present experiment. Sample 
temperature thus estimated within 68 sec [S3 after the 
start of bombardment was ~800°C for a beam power, 
P = 14 watt/cm 2 corresponding to the maximum flux 
(140/iA/cm 2 ) employed in the present case for bombard- 
ment at 100 keV. As the period of bombardment is much 
higher than this rise time, beam induced heating is ex- 
pected to be appreciable in the present case for the bom- 
bardment without scanning. Thus, thermally activated 
surface diffusion plays an important role here giving rise 
to a negative exponent I ~ E~°- 56 , as expected from lin- 
ear theory. It seems, the surface temperature profile on 
the sample generated by unscanned focused ion beam has 
made thermally activated diffusion process the dominat- 
ing one. This was not the case even when experiments 
were performed by increasing the bulk temperature of 
the sample |J. It should be mentioned here that we did 
not observe any flux dependence of the measured I in the 
flux range of 15 to 140 /iA/cm 2 used here as predicted in 
BH model. 

On the otherhand, for irradiation with beam scan- 
ning, effective beam intensity on the bombarded spot is 
reduced as the beam is sweeped continuously over the 
sample. The magnetic beam sweeping system of our im- 
planter uses scan frequencies of 0.51 Hz and 5.3 Hz along 
the horizontal and vertical direction to maintain beam 
homogeneity on the irradiated area of the sample. This 
means, each region of the bombarded spot of the size of 
beam diameter will be exposed to the ion beam every 
after few seconds only which is much less than the time 
of rise of the maximum equilibrium temperature due to 
beam heating. So promotion of surface diffusion of ther- 
mal origin will be hampered and ion induced surface dif- 
fusion will be the main relaxation process in the case 
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of beam scanning. However, at relatively lower energy 
thermal contribution may become noticeable. 

We can rewrite the ripple wavelength expressed by 
Eq. 2 by splitting the two diffusion terms as 



M 

which can be approximated as 

I 2 = A + BE- 1 +CE 2 



(4) 



(5) 



where A, B and C are fitting parameters. In this approx- 
imation, we have assumed that both depth (a) and lon- 
gitudinal spread (a) are linear with energy, as obtained 
exponents from TRIM calculations [H are close to unity 
(refer Fig. 2). It is possible to determine the shape of 
the distribution of the deposited energy density in the 
near surface region induced by the ion collision cascades. 
As for example, for the case of ripple topography for- 
mation on Ar bombarded graphite in 2 -50 keV energy 
region Habenicht et al Q obtained a power law relation 
1(E) ~ E p (p — 0.95) while a and a scale with ion en- 
ergy implying an energy-independent lateral spread of 
the damage cascade. However, we obtained 1(E) ~ E p 
(p = 0.45) with beam scanning (Fig. 2). A lower ex- 
ponent (p <C 1) is an indication of the changing lateral 
spread of the damage cascade with the ion energy for 
the energy region selected in the present experimental 
situation. The modified expression (Eq. 5) of the en- 
ergy dependence of ripple wavelength now considers the 
simultaneous effect of ion induced effective surface diffu- 
sion and thermally activated surface diffusion. 

The growth and decay of / with E of the present exper- 
imental results (Fig. 3a) can be fitted with Eq. 5 with fit- 
ting parameters A = 311697, B = 2.5 x 10 6 and C = 41.5 
for the ripple formation with beam scanning. The corre- 
sponding values for the ripple wavelength without beam 
scanning is A = -195533, B = 2.9 x 10 7 and C = 8.0. 
Here I and E are expressed in nm and keV respectively. 
The higher value of the coefficient of the second term for 
fitting the ripple wavelength data without beam scanning 
signifies that contribution of thermally activated surface 
diffusion is higher than that prevalent with beam scan- 
ning. 

To test the validity of our proposed formalism we have 
also used it for the data of Flamm et al jl8| where rip- 
ple structures with a characteristic wavelength between 
30 and 300 nm are formed on fused silica at room tem- 
perature after irradiation by broad beam ion (Ar) source 
in the energy region 0.6 to 1.5 keV. The expected domi- 
nance of ion induced effective surface diffusion has been 
reflected in this case by the excellent fit of the data 
(Fig. 3b) with Eq. 5 having much higher value of the 
coefficient of the ion induced diffusion term (C — 19861) 
over the value of the coefficient (B = 6077.8) represent- 
ing the thermal diffusion term. Thus, our proposed for- 



mula behaves consistently with homogeneous bombard- 
ment condition from the low (~ 1 keV) to the medium 
(~ few hundreds keV) energy regime where the concept 
of linear collision cascade model holds for heavy ion bom- 
bardment on solid target. 

In summary, we have experimentally verified, for the 
first time, the existence and interplay of two types of 
diffusion mechanism prevalent in ion induced nanoscale 
nanoscale ripple formation giving rise to the increasing 
or decreasing tendency of the ripple wavelength with ion 
energy. We have proposed a formalism for the ripple 
wavelength variation with the ion energy and verified its 
applicability over a wide energy range. The proposed 
formalism indicates that both thermal and ion induced 
surface diffusion act simultaneously to form stable ripple 
morphology. Ion beam induced surface temperature pro- 
file, rather than bulk temperature of the sample, seems 
to initiate the thermally activated diffusion process. It 
will be interesting to investigate the role of this tempera- 
ture profile in the formation of uniform nanoscale ripple 
structure and in the nonlinear phenomena. 

We would like to thank Mr. S. Roy and Mr. A. Das for 
their technical assistance and Dr. S. Hazra for valuable 
discussions. 
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FIG. 1: AFM micrographs of 40 Ar + bombarded Si surfaces 
at two representative energy, namely, 50 keV and 100 keV : 
(a) and (b) with ion beam scanning; (c) and (d) without beam 
scanning. Arrow marks indicate the projection of ion beam 
onto the surface. Variation of mean amplitude of the ripples 
with ion energy is shown in the lower plot. 



FIG. 2: Left scale depicts the variation of the mesaured wave- 
lengths I as a function of the 40 Ar + energy and the right scale 
shows the TRIM calculations for depth a and longitudinal 
straggling a as a function of ion energy. Power law fits to the 
wavelengths and TRIM data are inserted to each plot. The 
error bar of the measured data points have been determined 
after several repeated experiments and each data point rep- 
resents the averages over different regions of the bombarded 
spot. 



FIG. 3: (a) The square of the measured ripple wavelengths of 
the present experiment are fitted with Eq. (4). The function 
with the fitted parameters for the ripple wavelength data ob- 
tained with and without beam scanning are inserted within 
the graph, (b) The fitting of the data extracted from the 
measurements of Flamm et al [Q 
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